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The atomic dipoles in the Rabi lattice model exhibit quantum Ising dynamics which is exact in 
the limit of strong atom-photon interaction. It governs the para- to ferro-electric phase transition 
in the ground state. As a result of the strong coupling quantum Ising dynamics, the Rabi lattice 
model on an open chain also realizes two Majorana fermions localized at the boundary ends. In one- 
dimension, the DMRG calculations of the quantum phase diagram and atomic polarization clearly 
support the strong coupling behavior, and also yield the correct end-to-end dipole correlation which 
is a definite signature of these Majorana fermions. 
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Introduction. — In recent years, the lattice models of 
quantum cavities have received much attention for emu- 
lating the Bose-Hubbard model by exhibiting the Mott- 
insulator to superfluid type quantum phase transition for 
photons [1-5]. A quantum cavity (also called a cavity- 
QED system) refers to a strongly interacting matter- 
radiation system inside a high-Q resonator (cavity) . The 
simplest of it can be modeled as a single two-level atom 
interacting with a single mode of quantized radiation. 
It has been variously discussed that the lattice models 
of such quantum cavities may be realized, for example, 
by engraving an array of cavities in a photonic band- 
gap material [3, 6], or as a circuit-QED system [7-9]. 
It is hoped that such engineered systems will be useful 
in studying interesting quantum many-body problems, 
much like what the cold atoms in optical lattices do [10]. 

In complex physical systems, one is often concerned 
with the phase transitions in the ground state. For exam- 
ple, the superfluid to insulator transition in the ground 
state of the Bose-Hubbard (BH) model is such a phase 
transition driven by the competition between the kinetic 
energy and the local repulsion. Another very famous case 
which beautifully typifies quantum phase transitions is 
the quantum Ising (QI) model [11-13]. It describes an 
Ising system in an external field applied transverse to the 
Ising direction. The QI model has not only been used to 
study a variety of physical problems [11, 14-16], but it 
has also served as an important point of reference in the 
general understanding of quantum phase transitions [17]. 

The one-dimensional (Id) spin- 1/2 QI model with only 
nearest-neighbor interactions is an exactly solvable prob- 
lem that has been rigorously worked out by Pfeuty using 
the Jordan- Wigner fermionization of the spin-1/2 oper- 
ators (a method pioneered by Lieb, Schultz and Mattis 
for the XY chain) [13, 18]. For the Id QI model, in the 
fermionized form, Kitaev made a remarkable observation 
that the two (boundary) ends of an open chain carry a 
Majorana fermion each [19]. A fermion which is an anti- 
particle of itself is called a Majorana fermion [19-21]. 
While there has always been a great interest in finding 
the Majorana fermions in nature, the developments in 



quantum computation and condensed matter have fur- 
ther invigorated their search, as clear from the recent 
experimental activities [22]. Their proposed role in the 
intrinsically error-free (topological) quantum computa- 
tion indeed makes the realization of Majorana fermions 
a very valuable pursuit [19, 23]. 

In this Letter, we present an interesting case for the 
realization of Majorana fermions, through the QI dynam- 
ics, in the Rabi lattice model. By Rabi lattice we mean 
a lattice of the Rabi quantum cavities where each cav- 
ity has a two- level atom (or a spin-1/2) interacting via 
dipolar interaction, p ■ E, with a single mode of quan- 
tized radiation, and the inter-cavity coupling amounts to 
photon hopping. We show that the Rabi lattice model 
in the limit of strong atom-photon interaction rigorously 
tends to the QI model, which of course implies that the 
Rabi lattice in Id must have two Majorana fermions. 
This is in contrast to the commonly studied lattice mod- 
els formed by the Jaynes-Cummings type cavities, which 
show superfluid-insulator quantum phase transition for 
photons, but do not exhibit the QI behavior. Below we 
systematically discuss the emergence of QI dynamics in 
the Rabi lattice model, and support it by the density ma- 
trix renormalization group (DMRG) calculations in Id. 

Rahi lattice model. — A Rabi quantum cavity is a min- 
imal matter-radiation problem of great usefulness. It is 
described by the Hamiltonian: = u>{a'd + 5) + §cr^ + 
ja^{a^ + a), called the Rabi model. Here, uj is the pho- 
ton energy and e is the atomic transition energy. The 
dipole interaction, p-E, is written as 7(7^(0^ -I- a), where 
7 denotes the atom-photon coupling, d (a^) is the photon 
annihilation (creation) operator, and the Pauli operator, 
CT^ = \e){g\ + |5)(e|, measures the atomic dipole moment, 
p [24]. Moreover, = |e)(e| — \g){g\. The kets \g) and 
|e) denote the atomic ground and excited levels, respec- 
tively. The simplicity of Hfi is deceptive as its eigen- 
value problem has not been solved exactly, although it 
has been studied extensively [25-27]. A popular variant 
of this problem, due to Jaynes and Cummings, is however 
exactly solvable [24, 28] . 

In the Jaynes-Cummings (JC) model, the atom-photon 
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interaction is taken as 'y{a~^a + (T~a^), by dropping 
(T+a^ + a~d from the dipole interaction. This is called 
the rotating- wave approximation (RWA) . It is based on a 
dynamical argument that, close to resonance (w ~ e), the 
processes corresponding to ct^o^ + a~a are faster com- 
pared to (T+a -|- a^a\ and hence can be neglected. Here, 
(T+ = \e){g\ and = The RWA, however, dras- 

tically alters the symmetry properties which gives rise 
to qualitatively different physical behavior in the two 
models. In the JC model, it leads to the conservation 
of 'polariton' number, a^a + a^a~ , and makes it ana- 
lytically solvable. This number conservation [continuous 
U{1) gauge-symmetry] is not present in the Rabi model. 
Instead, it is discretely symmetric under the parity op- 
erator, — xcr^, where x — — e*''^" and n — a)a. 
The parity is -(-1 for even polariton number and —1 for 
the odd. Thus, in a Rabi cavity, the states with even 
polariton number do not mix with the odd ones. 

In the lattice models of such quantum cavities, the 
photons can coherently hop from one cavity to another, 
while each atom stays put inside the cavity it belongs to. 
This hopping is a process different from the incoherent 
decay of photons in a cavity. The photon does not decay 
while hopping, instead it disperses over the lattice. The 
time scale of incoherent losses is assumed to be much 
longer than the time scales associated with the photon 
hopping and the atom-photon coupling, 7. 

Notably, the studies that find the Bose-Hubbard type 
quantum phase diagram for photons have primarily fo- 
cussed on the Jaynes-Cummings lattice models. As noted 
above, the RWA introduces polariton number conserva- 
tion, due to which the JC lattice model exhibits polari- 
ton condensation, just like the material particles of the 
Bose type. It would be interesting, therefore, to inves- 
tigate as to what happens without the RWA. Inside a 
high finesse cavity with strong atom-photon interaction, 
the counter-rotating terms would invariably influence the 
physical properties. After all, they are alway there in the 
dipole interaction. Hence, we study the lattice model of 
the Rabi quantum cavities. 

The Rabi lattice model can be written as follows. 
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Here, Hrj — uj{fii -t- 5) -t- |crf -t-7(Tf (a| -|-a;), is the 'local' 
Hamiltonian of the Z*'' Rabi cavity in the lattice. The 
second term in H describes the nearest-neighbor photon 
hopping with (amplitude t), and 5 denotes the nearest- 
neighbors (without counting any nearest-neighbor bond 
twice). For example, on square lattice, 5 — ax and ay 
(and not ±ai and iay). One can also include further 
neighbor hopping. But, we only consider the nearest- 
neighbor case. Moreover, the parameters u, e, 7 and t 
are taken to be positive real numbers. 

Note that, for the Rabi lattice model [H of Eq. (1)], the 
parity of the individual cavities is not conserved anymore. 



This is because the photon hopping incessantly changes 
the photon occupancy, fii, in each cavity which causes 
the the local parity to fluctuate between odd and even. 
However, the global parity defined as, (— )^ Hi Xi'^'f ' 
still a symmetry of the H. Here, L is the total number 
of lattice sites (cavities). Since the Rabi lattice model 
only has a global parity (discrete) symmetry, a thermo- 
dynamic transition, if it occurs, would be characterized 
by the breaking of it, and would not be accompanied by 
gapless excitations in the 'ordered' phase. This point has 
been noted recently for the Rabi lattice model [29]. In 
fact, we expect the Rabi lattice model to exhibit quantum 
Ising transition, as hinted at by the fluctuating local par- 
ity picture presented above. In the following, we clearly 
establish for this physical scenario by suitably transform- 
ing the H to & form in which the quantum Ising dynamics 
expresses more evidently, and explicitly corroborate it by 
the DMRG calculations in Id. 

Quantum Ising dynamics. — Consider the atom-photon 
interaction, erf (a|-|-a;), in a Rabi cavity. Since the atomic 
dipole operator af has eigenvalues ±1, we can absorb it 
into the electric field, (a| + &;), by the unitary transfor- 
mation under U = H^^i where Ui ~ P^^ + P; x; and 
= (1 ± erf )/2. Under this U, we get, U'^aiU = afai, 
U'fniU = hi, U^a^U = x/crf , and U^afU = af. The op- 
erators hi and af are obviously invariant under U. More- 
over, the a; and aj carry af, as wanted. Due to this, the 
dipole interaction transforms to a static 'displacement' 
term, j{aj -I- O;). Furthermore, the local parity, — x;crf, 
transforms to —erf. Thus, in the transformed version of 
the Rabi cavity, the dipole interaction turns into a dis- 
placement field for photons, the erf is to be understood 
as parity, and erf continues to be the atomic dipole. 

Now consider the Rabi lattice model, H. Under U, it 
transforms to the following form. 
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The role of hopping in causing local parity fluctuations 
is more evident now due to the explicit (^fcrf^g factors. 
Although the transformed Hamiltonian, Hi, appears to 
be more complex, it actually helps in making systematic 
progress. Look at the atom-photon interaction. It acts 
as a constant displacement field for the photons, which 
guarantees that (ai) ^ in the ground state. Thus, a 
static 'electric' field, proportional to 7, is ever present. It 
has interesting consequences for the parity and dipole dy- 
namics. This static field, through hopping, would gener- 
ate a static 'Ising' interaction between the atomic dipoles, 
erf. Moreover, the atomic energy, e (renormalized by the 
expectation of x/), would act as a 'transverse' field on the 
parity, erf . Thus emerges the quantum Ising dynamics in 
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the Rabi lattice model. While the photon fluctuations 
above the static field would also be present, in the strong 
coupling limit (7 ^ the other energy scales), they can be 
neglected. 

To give this discussion a proper theoretical form, we do 
a displacement operation on Hi such that the 7(0} + ai) 
is absorbed into hi. This is done by the unitary operator, 

the Hi transforms 
to the following form under D. 



D = nr=i A, where Di 
ai- ^ and D^xiD = X/ e 
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This is the Rabi lattice model written explicitly in terms 
of the static field, J, and the photon fiuctuations. Now, 
invoke the strong coupling limit, ^ S> 1. It amounts to 
neglecting the photon operator terms in H2. The limiting 
strong coupling problem that results is the quantum Ising 
model of parities and dipoles written below. 
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Here, the Ising interaction, J = —2t^, between atomic 
dipoles is ferro-electric, and the transverse field acting 
on the local parities, h = |exp(— -^), is exponentially 
suppressed (akin to the Ham, or Frank-Condon, reduc- 
tion factors in vibronic states [30, 31]). 

Ground state properties. — Since the QI model is a well- 
studied problem, having Hqj as the limiting Hamiltonian 
immensely helps in understanding the ground state of the 
Rabi lattice model. This is more so when 7 is the largest 
energy scale (the domain of applicability of Hqj), but 
also in general. It predicts two distinct phases, character- 
ized by the atomic polarization, that undergo quantum 
Ising transition between them. 

For I J| > h, the Rabi lattice would exhibit spontaneous 
polarization. In this ferro-electric (FE) ground state, the 
order parameter, p = ^li'^i) 0. For | J| < /i, each 
Rabi cavity on the lattice would (roughly) behave as in- 
dependent, and p = [para-electric (PE) phase]. This 
can also be restated in terms of photons, as often done in 
the literature on the cavity lattice models. The photon 
order parameter, ip — (a/), can be described (after hav- 
ing applied U and ID) asijj = {(jf{ai — ^)) ~ ~Pz]- Thus, 
in the FE phase, the photons exhibit 'superfluidity'. The 
PE phase may likewise be called a Mott insulating phase. 
However, it not quite the same because a Mott phase is 
characterized by an integer polariton number per cavity. 



which is not true for the Rabi lattice. Moreover, the FE 
phase is also not the usual superfluid described by a com- 
plex U{1) order parameter with gapless modes. Instead, 
it has a real Ising (Z2) order parameter with gapped ex- 
citations. Notably the excitations in both the phases are 
gapped. They can be viewed as parity flips in the PE 
phase, and atomic dipole flips in the FE phase. The clos- 
ing of the gap marks the quantum phase transition. 

The factor, exp(— -^^a-), in the transverse field, h, too 
has some important consequences. For ^ 3> 1, it renders 
the variations in e unimportant, as h would invariably be 
small regardless of the value of e. For the same reasons, 
it is immaterial whether e is near resonance {'^ uj) or far. 
The smallness of h also allows the Ising interaction, J, 
to dominate because of which the critical hopping, tc, for 
the transition from PE to FE phase would be small. 

To explicitly validate the strong coupling behavior, and 
to see deviations from it, we numerically investigate the 
Rabi lattice model in Id using DMRG. The DMRG is a 
very nice method for studying the ground state proper- 
ties of Id quantum problems [32, 33]. Since the Hqj in 
Id is exactly solvable, it also allows us to make proper 
comparisons with the DMRG results. Notably, through 
the strong coupling QI problem, the Id Rabi lattice also 
presents an exciting possibility for realizing Majorana 
fermions, which indeed makes it a special case to study. 

For Hqj in Id, the exact polarization in the FE phase 
is p = {l~ js)^ (in the thermodynamic limit) [13], where 

X = \ J\/h ^ f^exp{^). Using DMRG, we compute 
p as a function of t in the ground state of the Rabi lattice 
model (in the H2 form). It clearly exhibits the PE to 
FE phase transition, and agrees nicely with the strong 
coupling QI form, as shown in Fig. 1. This DMRG data 
is computed for a chain of length L — 204, and by taking 
four photon states (n; = 0, 1, 2, 3) per Rabi cavity. 

Since the exact quantum critical point for the PE to 
FE phase transition in Hqj in Id is A = 1 [13], the strong 

coupling critical hopping is tc = |^ exp (— 735-)- The tc 
for e = UJ that has been derived recently in Ref. 34 is just 
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FIG. 1. Polarization, p, vs. hopping, t, in Id. The black 
continuous line is the exact strong coupling prediction. 
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FIG. 2. Quantum phase diagram for the Rabi lattice model 
in Id. The PE (FE) denotes the para- (ferro-) electric phase. 
The continuous line is the strong coupling phase boundary. 

a particular case of the present result. Note that, for a 
given ^, the strong coupling tf. scales as e, and it is true 
in any dimension. This implies that, in the plane of j 
and ^, the strong coupling phase boundary would be a 
universal curve, independent of e. In Fig. 2, the DMRG 
calculated phase boundary is compared with the strong 
coupling prediction. Clearly, for — ^ 1, the numerical 
data approaches the universal tc/e derived above. 

For ^ < 1, however, there are clear deviations from the 
strong coupling QI behavior. We can approximately ac- 
count for these deviations within a simple-minded impro- 
visation of the Hqi , which corresponds to replacing oj by 
uj~ztp\ in Eq. (4). It is basically correcting w for the pho- 
ton dispersion. Here, z is the nearest-neighbor coordina- 
tion (which is 2 in Id), and p\ = {(jfaf_^_^) is the nearest- 
neighbor correlation. This improvisation is effected by a 
more general displacement operation, Y\i^^°''"°''-^°''\ on 
Hi, and approximating ta^i ^ afaf_^g{ai + ai-fs + h.c.) 
(arising from the photon hopping after applying the dis- 
placement) by ztapf + ^')- Putting the terms lin- 
ear in photon operators to zero gives a = j/{u — ztpf), 
and ignoring the photon fluctuations gives the improvised 
Hqj. Thus, we get the improvised QI phase boundaries 
shown in Fig. 2, whose agreement with the DMRG data 
is quite reasonable. Approximating erf erf by is par- 
ticularly okay near the transition where pi is strong (sup- 
ported by DMRG calculations of pf). At the transition 
point, Pi — ^ for the Id QI model. The improvised QI 
model may not work for much weaker couplings. 

Majorana fermions. — On an open chain, the strong 
coupling Hqi guarantees that the Rabi lattice has two 
'unpaired' Majorana fermions in the FE phase. This is 
because the exact excitations of the QI model in Id be- 
have as fermions, of which, one corresponds to having 
two Majorana fermions localized at the two ends in the 
ordered phase [13, 19]. Another important feature con- 
cerns the exact end-to-end correlation, p^j^ = (crfaf^) = 

+ 0{1/L), in the ground state of the open chain QI 
model [13]. The fact that pfj^ is (in the thermody- 
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FIG. 3. End-to-end correlation, pj^, vs. on an open chain. 



namic limit) is at odds with the normal expectation, p^. 
This striking property of pfj^ is a direct consequence of 
these Majorana fermions [18], and thus a definite signa- 
ture thereof. Our calculations clearly find this behavior 
in the Rabi lattice model for large couplings (see Fig. 3). 
From an experimental perspective. Fig. 3 is quite inter- 
esting because it compares two observables only (end- 
to-end dipole correlation and spontaneous polarization), 
and doesn't explicitly involve any model parameters. 

Summary. — The properties of the Rabi lattice model 
crucially depend upon the dipole interaction and photon 
hopping. For the dipolar atom-photon interaction, the 
parity is the only good quantum number locally in each 
cavity. The photon hopping, however, causes these local 
parities to fluctuate. This is aptly described by Hi, in 
which the atomic dipoles, through photon hopping, com- 
pete with the parities. This identification is facilitated 
by the unitary transformation, U, which helps the pho- 
tons to carry the discrete Z2 phase of the dipoles. It 
also transforms the atom-photon interaction to a static 
displacement field for photons. This leads to a crucial 
insight that if the atom-photon coupling is very strong, 
then the photons will behave as a constant 'electric' field 
with negligible fluctuations. This is effected by rewriting 
the problem in the displaced photon basis as H2, which 
in the strong coupling limit, ^ » 1, tends to Hqi, the 
quantum Ising model of dipoles and parities. This drives 
the PE to FE quantum phase transition, characterized 
by the spontaneous polarization, p, of atomic dipoles, in 
the Rabi lattice model. In the FE phase on an open 
chain, the strong coupling QI dynamics guarantees the 
existence of two Majorana fermions located at the two 
ends. This is an important finding in view of the current 
excitement about Majorana fermions. The end-to-end 
dipole correlation equals p^ is a clear signature of the 
Majorana fermions to look for experimentally. The theo- 
retical analysis of the Rabi lattice model is well supported 
by the numerical DMRG calculations in Id. 
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